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ABSTRACT 

The supergiant fast X-ray transient IGR J08408-4503 was discovered by INTEGRAL on May 15, 
2006, during a bright flare. The source shows sporadic recurrent short bright flares, reaching a peak 
luminosity of 10 36 erg s~ x within less than one hour. The companion star is HD 74194, an Ob5Ib(f) 
supergiant star located at 3 kpc in the Vela region. We report the light curves and broad-band spectra 
(0.1-200 keV) of all the three flares of IGR J08408-4503 detected up to now based on INTEGRAL 
and Swift data. The flare spectra are well described by a power-law model with a high energy cut-off 
at ~ 15 keV. The absorption column density during the flares was found to be ~ 10 21 cm -2 , indicating 
a very low matter density around the compact object. Using the supergiant donor star parameters, 
the wind accretion conditions imply an orbital period of the order of one year, a spin period of the 
order of hours and a magnetic field of the order of 10 13 G. 

Subject headings: X-rays: binaries - X-rays: bursts - X-rays: individual: IGR J08408-4503 



1. INTRODUCTION 

Supergiant Fast X-ray Transients (SFXTs), consisting 
of a wind accreting compact object and an OB super- 
giant donor star, spend most of the time in a quiescent 
state, with X-ray luminosities of the order of 10 32 — 10 33 
erg s _1 . Sporadically they go into outburst reaching lu- 
minosities of 10 36 — 10 37 erg s^ 1 (Sguera et al. 2005; in't 
Zand 2005; Negueruela et al. 2006; Smith et al. 2006). 
The outbursts or strong flaring activities are character- 
ized by very short time scales, i.e. from minutes to hours. 
SFXTs have been discovered only recently as a new class 
of X-ray sources, and may unresolved fundamental ques- 
tions remain open. The lack of confirmed orbital or pulse 
periods complicates the determination of the system pa- 
rameters as well as of the compact object nature. Up 
to now IGR J11215-5952 alone shows regularly recur- 
rent outbursts every ~ 330 days, most likely linked to 
its orbital period (Sidoli et al. 2006). To account for 
the accretion process that causes the short outbursts, it 
was proposed that the clumpy wind in early-type stars 
could be captured by the compact object producing the 
X-ray flares on the observed time scale (in't Zand 2005). 
However, this scenario has yet to be well understood and 
confirmed for OB stars in binary systems. 

The number of known SFXTs has grown recently 
thanks to the International Gamma-Ray Astrophysics 
Laboratory {INTEGRAL), see e.g. IGR J16465-4507, 
IGR J17544-2619, or IGR J11215-5952 (e.g. Walter et 
al. 2006; Sguera et al. 2006; Sidoli et al. 2006). The new 
SFXT presented here, IGR J08408-4503, was discovered 
in the Vela region on May 15, 2006 with INTEGRAL 
during a bright outburst lasting about 15 minutes (Gotz 
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et al. 2006). The candidate optical counterpart was first 
tentatively identified as the supergiant, Ob5Ib(f), HD 
74194 star (Gotz et al. 2006; Masetti et al. 2006) lo- 
cated at 3 kpc in the Vela region (e.g. Walborn 1973; 
Humphreys 1978; Schroder et al. 2004). Successively, us- 
ing the Swift X-ray Telescope (XRT) a refined source po- 
sition (ajaooo = 08 h 49 m 47 s .97 and 5 J2 ooo = -45 o 03'29'.'8 
with an uncertainty of 5"4) was derived, strengthening 
the association with HD 74194 (Kennea & Campana 
2006). 

Optical spectra of HD 74194 measured a few days af- 
ter the outburst (Barba et al. 2006) strongly resemble 
the ones of the companion star of another SFXT, IGR 
J17544-2619 (Pellizza et al. 2006). In addition the ra- 
dial velocity variations in the Hel and Hell absorption 
lines with an amplitude of about 35 km s _1 were mea- 
sured, suggesting either a pulsating variable nature of the 
supergiant star HD 74194 or a possible Doppler orbital 
modulation (e.g. Conti et al. 1977). 

IGR J08408-4503 shows the typical SFXTs recurrent 
short flaring events. Using INTEGRAL archival data, 
an outburst was detected on July 1, 2003 (Mereghetti 
et al. 2006), while the last bright flare was observed on 
October 4, 2006 with the Swift satellite (Ziaeepour et 
al. 2006). 

In this letter we present the light curve and broad-band 
spectral study of all the three outbursts of IGR J08408- 
4503 observed up to now. In addition, thanks to the 
knowledge of the companion star parameters, the possi- 
ble wind accretion conditions are discussed. 

2. DATA ANALYSIS AND RESULTS 

IGR J08408-4503 was observed twice in outburst with 
INTEGRAL (Winkler et al. 2003), on July 1, 2003 and 
May 15, 2006, i.e. during the satellite orbits 83 and 438. 
This dataset includes publicly available data and part of 
A03 Vela region guest observations. Hereafter we call 
the two datasets flare 1 and flare 2, respectively. We 
used the data from the coded mask imaging telescope 
IBIS/ISGRI (Ubertini et al. 2003; Lebrun et al. 2003) 
at energies between 15 keV and 200 keV and from the 
JEM-X monitor (Lund et al. 2003) between 3 and 20 
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keV. For JEM-X the data were extracted only for flare 
2, flare 1 being outside the JEM-X field of view. The 
data reduction was performed using the standard Offline 
Science Analysis (OSA) version 5.1. 

Single pointings (~ 2000 s each) were deconvolved and 
analyzed separately. The IBIS/ISGRI 15-40 keV high 
energy light curve, shown in Fig. 1, has been extracted 
from the images using all available pointings. Light 
curves with a 100 s time bin were also extracted around 
the outburst peaks. 

The third outburst, flare 3, was observed on October 4, 
2006 (14:45:42 UT) with Swift/BAT (Gehrels et al. 2004; 
Barthclmy ct al. 2005). Since the flare was detected with 
an image trigger (15-150 keV), only "Survey" data prod- 
ucts are available from the BAT instrument, with a typ- 
ical integration time of ~ 300 s. We note that the BAT 
data coverage is not continuous due to the low Earth 
orbit of the Swift satellite. For each ~ 300 s time bin 
the target count rate in the 14-40 keV energy range was 
evaluated with the mask-weighting technique using the 
XRT coordinates. The BAT data were selected 15 h be- 
fore and 30 h after the trigger, and reduced using the 
standard software version 2.5. 

Swift slewed automatically to the direction of the 
source and XRT (Burrows et al. 2005) (0.1-10 keV) 
follow-up observations were performed. The XRT data 
were reduced using ftools version 6.1.1., and processed 
with XRTPIPELINE version 0.10.4, using standard filter- 
ing selections. The XRT 0.1-10 keV energy band light 
curves have been extracted by selecting a circular region 
of 20 pixels around the source for the Photon Count- 
ing (PC) mode, and rectangular region of 40x20 pixels 
for the Window Timing (WT) mode. Background light 
curves were extracted in source free regions in order to 
produce the final background-subtracted XRT/PC/WT 
light curves. The count rate of the source in both modes 
is low enough in order to avoid the pile- up in the CCD. 

2.1. Flare light curves 

In Fig. 1 we show the three outburst light curves 
of IGR J08408-4503 observed up to now. The light 
curves are dominated by a first bright flare that typi- 
cally reaches a peak luminosity of 3 — 6 x 10 36 erg s" 1 
(see Sec. 2.2) in a time scale of several tens of min- 
utes to hours. The observed flares vary from an isolated 
strong flare (flare 2) to a more structured flaring activ- 
ity (flare 1 and 3). Moreover flare 3 shows also weaker 
flares with a luminosity of the order of 5 x 10 35 erg s _1 . 
The three outburst durations are of the order of sev- 
eral hours. The flare characteristics are similar to those 
observed in other SFXTs, such as XTE J1739-302 (e.g. 
Sguera et al. 2005), IGR J16465-4507 (e.g. Lutovinov ct 
al. 2005), IGR J17544-2691 (e.g. in't Zand 2005), and 
IGR J11215-5952 (e.g. Sidoli et al. 2006), but resemble 
also those well studied in persistent wind accretion high 
mass X-ray binary (HMXB) sources like Vela X-l (e.g. 
White et al. 1983), 4U 1700-377 (e.g. Rubin et al. 1996). 

Up to now only IGR Jl 1215-5952 has shown a - 330 
days periodic outburst. Therefore, using a linear orbital 
function T(n) = Tq + nP or h we attempted to test the 
hypothesis of a periodically recurrent outburst also for 
IGR J08408-4503. P orb is the orbital period in days and 
T is taken from flare 2 (MJD 53780.771215). A possible 
orbital period is 4.2988 days (with n\ = 244 and n 2 = 



33), however, searching in the INTEGRAL archival data 
around the expected outburst time, we did not find any 
significant flux excess at the IGR J08408-4503 source 
position. We conclude that the three observed outbursts 
are sporadic episodes of accretion of matter fed from the 
wind of the supergiant companion star. 

In order to directly compare the light curves measured 
with different instruments (see Fig. 1), all the count rates 
have been converted to the 0.1-10 keV fluxes assuming 
the spectral shape derived in Sec. 2.2. 

2.2. Flare spectra 

We performed the spectral analysis using XSPEC ver- 
sion 11.3, for the ISGRI data (20-200 keV), flare 1, for 
the 3.5-20 keV JEM-X data with the simultaneous 20- 
200 keV ISGRI data, flare 2, and for the 0.3-100 keV si- 
multaneous XRT/BAT data, flare 3. Only the brightest 
parts of the flares have been considered. For the simul- 
taneous data a constant factor was included in the fit to 
take into account the uncertainty in the cross-calibration 
of the instruments. All data were rebinned in order to 
have 3cr points and the spectral uncertainties in the re- 
sults are given at a 90% confidence level for a single pa- 
rameter. We use the source distance of 3 kpc throughout 
the paper. 

The XRT/BAT (flare 3) broad-band (0.3-100 keV) 
dataset provides us the widest energy coverage and the 
best statistics, so we fitted it first. We fitted the 
data using a simple photoelectrically-absorbed power- 
law (pl) model which was found inadequate with a 
X 2 /d-o.f. = 132.76/76. The addition of a high-energy 
exponential cut-off significantly improved the fit to 
X 2 /d.o.f.=80. 36/75, resulting in a best-fit photon index 
of 0.1±0.2 and a cut-off energy at 15±5 keV. The hydro- 
gen column density, Ah, was found to be 1.0 ±0.3 x 10 21 
cm~ 2 . This value is compatible with Galactic value in 
the direction of the source, Ah = 3 x 10 21 cm~ 2 , re- 
ported in the radio maps of Dickey & Lockmann (1990). 
We attempted to fit the low energy spectrum also with 
a thermal black-body model, BB, plus a PL for the high 
energies, but we found high bb temperature, fcTsB ~ 8 
keV, which is much larger than the measured values for a 
neutron star (NS) surface or polar cap thermal emission. 
Also a fit with a thermal bremsstrahlung model can be 
statistically ruled out by our data (x 2 /d.o.f.=183.2/76). 
The energy spectra during the flaring activity are hence 
best fitted with a high-energy exponential cut-off PL 
model similar to the spectra observed from persistent 
wind accreting HMXB hosting a NS, e.g. Vela X-l 
(White et al. 1983). 

The joint JEM-X/ISGRI (3.5-200 keV), flare 2, spec- 
trum was also fitted with the cut-off PL model used for 
the flare 3. The best fit values are similar to the ones 
of flare 3, see Table 1. Given that we were not able to 
constrain the Nh value (as the JEM-X bandpass starts 
above 3 keV) we fixed it to the value found from flare 3. 

For flare 1 only ISGRI data (20-200 keV) are available. 
The data can be well fitted with a simple PL model with 
a photon index, T, of 2.5±0.5. Note that fitting flare 2 
(ISGRI) and 3 (BAT) spectra in the same energy range 
gives consistent photon index values, namely T=2.8±0.3 
and F=2.3±0.4, respectively. The apparently soft spec- 
tra derived using only the data above 20 keV show how 
important it is, to have a broad energy coverage in order 
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TABLE 1 

IGR J08408-4503 spectral parameters using a power-law 

MODEL WITH AN EXPONENTIAL HIGH ENERGY CUT-OFF 



Fig. 1. — Bottom, from left: IBIS/ISGRI (flare 1,2) and Swift (flare 3) light curves. IBIS data are binned at ~ 2 ks. In the last panel 
the crosses represent the extrapolated BAT flux, while the diamonds and triangles the XRT/WT and PC data. The time axis is expressed 
in hours after T =52821.0, Ti =53870.0, and T 3 =54012.615 MJD, respectively. The Swift light curve shows also weak flares e.g. at t~3.5 
h IGR J08408-4503 is detected in BAT at ~7 o~ level. Top: the brightest part of each outburst is shown on a smaller time scale (100 s 
time bin for ISGRI). 

to better characterize the spectrum of these sources. To 
be coherent, we applied the same spectral model derived 
from the broad-band fit of flare 2 and 3, also for flare 
1, even if we had a smaller energy range. We fixed the 
slope of the power law index to the values found (see Ta- 
ble 1) and found the cut-off energy compatible with the 
one found for flares 2 and 3. 

In Fig. 1, last panel, weak flares are visible after the 
main peak. To infer the flux at very low mass accretion 
rate we extracted a XRT/PC spectrum in the 0.1-6 keV 
energy band at ~2.2 hr after the main flare, see Fig. 1. 
This is the last time interval where the source is clearly 
detected in the XRT images. We call this weak flare, flare 
4. Using the same model found for flare 3, the source 
was at a bolomctric luminosity of 5 x 10 35 erg s _1 , with 
a marginal indication of an increase of the absorption 
column density. 

The best fit parameters for the different flares are re- 
ported in Table 1. In Fig. 2, we show the broad-band 
unfolded spectrum and the residuals of the data to the 
PL with a high energy cut-off model for flare 3. 
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Fig. 2. — Combined XRT (crosses) and BAT (diamonds) spec- 
tra of flare 3. The best fit PL high energy cut-off model and the 
residuals (lower panel) are shown. 



3. DISCUSSION AND CONCLUSIONS 

IGR J08408-4503 is a recurrent SFXT with most likely 
sporadic episodes of accretion of matter from the wind 
of the Ob5Ib(f), HD 74194, supcrgiant companion star. 
The broad-band spectrum (0.3-200 keV) allowed us to 
perform an improved spectral analysis for the new source 
IGR J08408-4503 using XRT/BAT and JEM-X/ISGRI 
data. The best fit to the data required a PL with a 
high energy cut-off, see Table 1. The spectrum is typical 
for HMXBs hosting a NS (White et al. 1983) where the 
emergent radiation is presumably produced by accretion 
columns at the magnetic poles of the NS. The bolometric 
flare peak luminosity is of the order of 3 — 6 x 10 36 erg 
s , which are typical values for wind accreting HMXB. 
At low mass accretion rates, the source is observed at a 
bolometric luminosity of 5 x 10 35 erg s _1 . The source 
quiescent luminosity is in the order of ~ 2 x 10 32 erg s _1 
(Kcnnea & Campana 2006). 

During the outbursts the derived absorption of ~ 1 x 
10 21 cm -2 is compatible with the total Galactic absorp- 
tion in that direction as estimated from the HI maps, 
indicating that the source is not surrounded by large 
amounts of material. Some information about the distri- 
bution of matter around an X-ray binary can be inferred 
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TABLE 2 

Stellar parameters for 
hd 74194 derived from 
(Lamers et al. 1995) 



Parameter 


Value 


T eff (k) 


33000 


R* (Rq) 


14.5 


M* (M ) 


28 


log(L*/L ) 


5.3 


M„ a (Mq yr- 1 ) 


2.3xl0~ 7 


Voo (km s 1 ) 


2000±300 


u csc (km s _1 ) 


776 



„ , .. a Calculatcd using .Eq. 2. ... , ., , , „ 

from observations oi trie variations with orbital phase of 
its spectrum caused by absorption along the line of sight 
to the X-ray star. For IGR J08408-4503, the orbital pe- 
riod is not known, therefore using a simply spherically 
symmetric stellar wind model, we approximate the wind 
density (e.g. see Eq. (7) in Heap & Corcoran 1992) at 
the NS surface (a » R*, a is the binary separation) for 
a circular orbit: 



N H « 1 x 10 21 x 

My, 

10- 7 Moyr- 1 



R* 



U.5R 



2000 km s" 1 



cm 



(1) 

-2 



where M w is the HD 74194 mass loss rate and v x is 
the terminal wind velocity. The used stellar parameters 
for HD 74194 are reported in Table 2. We found the 
estimated Ah compatible with the derived absorption 
from the spectral fit. On the other hand, locating the 
companion star at a — 2_R*, and using a maximal wind 
acceleration parameter (3 = 1, Nh will change by a factor 
ln(i?*) « 27, i.e. Ah w 2.7 x 10 22 cm" 2 . The observed 
Ah hints towards a large binary separation and therefore 
a large orbital period. The measured absorption is in 
variance with what has been observed for other SFXT 
(see e.g. Walter et al. 2006). 

The companion mass loss rate, M w , was calculated for 
HD 74194 using the average relation for the galactic O 
or B supergiants stars (Lamers & Cassinelli 1999), 

log(M„) = -1-37 + 2.07 log(L»/10 6 )-l°g("°cfl2- 5 )- (2) 

The M w value is a factor 10-100 lower that what found 
for the persistent wind accretion HMXB, like Vela X-l or 
4U 1700-37. The free parameter to obtain a higher Nn 
value for IGR J08408-4503 is to have a higher mass loss 
rates, M w . Using Eqs. (2) and (2) we can see that if the 
terminal wind velocity is a factor 5 lower (as observed for 
an intrinsically obscured high energy source IGR J 163 18- 
4848 (Filliatre & Chaty 2004)), we would have a higher 
absorption of ~ 10 23 cm -2 , but such a low wind velocity 
is not observed for HD 74194. This indicates that during 
the accretion of matter for IGR J08408-4503 there is no 
accretion wake of dense matter surrounding the compact 
object (Blondin 1994). 

We will assume that all the material within the capture 
radius 

Ra.cc. = ~ ~ (3) 



vi + K 



',9.2 x 10 



8 / M acc 



1/2 



/ M* + M x 
V 28M Q 



1/3 



PJ h cm 



is accreted by the compact object (Bondi & Hoylc 
1944). M x is the mass of the NS, v x its orbital velocity, 
and M acc = L x c~' 1 r)~ 1 the mass accretion rate, where 
r\ ~ 0.2 is the accretion efficiency for a NS. The rate of 
mass captured is then given by M acc = M w R\ cc j (4a 2 ), 
with i? acc = 9.3 x 10 9 cm. Wc consider that v x >> v x 
and a » From the measured persistent emission 
Li~ 2 x 10 32 erg s -1 , we derive from Eq. (4) an orbital 
period of ~1.5 years. 

The condition for accretion to take place is that the NS 
magnetosphere radius is within the mass capture radius 
and the corotation radius, i.e. RMag < -Race and Ruag < 
R COI (Illarinov, A. & Sunyaev, 1974). If we set i? acc = 



Rear we have P? 



spin 



2.36 x 10 27 (1.4M X )^- 



7000 s. 



The magnetosphere radius is given by 



R Mag = OA^MWu-WMWpW, 



(4) 



where \x = BR"^ S . Setting the accretion condition 
Ruag — Rear = Race the magnetic field has to be 
B ~ 1.1 x 10 13 (P orb /lyr)- 2/3 G, for R NS = 10 6 cm. 
Assuming the orbital period derived above, the NS mag- 
netic field has to be of the order of 10 13 G. These are 
typical magnetic field values for young HMXBs hosting 
a NS, like e.g. Vela X-l (La Barbera et al. 2003). 

The low A^h value measured during the flares is not 
consistent with the picture in which they are caused by 
clumps in the donor wind. In an alternative scenario, 
the flares could be associated with the sudden accretion 
onto the magnetic poles of matter previously stored in 
the magnetosphere during the quiescent phase. However, 
in order to have such a mass storage, the above simplest 
accretion conditions have to be studied for different sce- 
narios (e.g. R acc > R COI > i? mag ). One can derive these 
conditions by varying opportunely P rb, -Pspin, and B 
(Bozzo et al. 2007). 

We conclude that these recurrent sporadic very short 
outburst episodes, due to the accretion of matter from 
the wind of a supergiant companion star, imply a spin 
period of the order of hours with a long orbital period, 
and a 10 13 G magnetic field for the NS. The determina- 
tion of all the system parameters can help to solve the 
accretion mechanism. 
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